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Abstract. 

We present a general method for analytically factorizing the n-fold form factor integrals 
fnlrV) f° r me correlation functions of the Ising model on the diagonal in terms of the 
hypergeometric functions 2F1 ([1/2, N + 1/2]; [iV + l];i) which appear in the form factor 
ftiJi W' New quadratic recursion and quartic identities are obtained for the form factors for n = 
2, 3. For n = 2, 3, 4 explicit results are given for the form factors. These factorizations are 
proved for all N for n = 2, 3. These results yield the emergence of palindromic polynomials 
canonically associated with elliptic curves. As a consequence, understanding the form 
factors amounts to describing and understanding an infinite set of palindromic polynomials, 
canonically associated with elliptic curves. From an analytical viewpoint the relation of 
these palindromic polynomials with hypergeometric functions associated with elliptic curves 
is made very explicitly, and from a differential algebra viewpoint this corresponds to the 
emergence of direct sums of differential operators homomorphic to symmetric powers of a 
second order operator associated with elliptic curve. 
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1. Introduction 

The form factor expansion of Ising model correlation functions is essential for the study of 
the long distance behavior and the scaling limit of the model. This study was initiated in 
1966 when Wu [1] computed the first term in the expansion of the row correlations both for 
T > T c , where the result is a one dimensional integral, and for T < T c , where the result is a 2 
dimensional integral. By at least 1973 it was recognized [2 1 that the diagonal correlations and 
form factors are a specialization of the results for the row correlations. The extension to form 
factors for correlations in a general position and from the leading term to all terms was first 
made in 1976 [3|. This leads to the general result that for the two dimensional Ising model 
with interaction energy £ = - Yjj,k{E v °~ j^j+i ,k + E h cr j yk cr j^+i), with cr^j = +1, the form 
factor expansion for T < T c is 



(cr Q , Q cr M , N ) = (1 - f) 1/4 ■ {1 + J] O' (1) 
where? = (sinh2£ v /l i; :r sinh 2E h /k B Ty 2 , and for T > T c 

00 

{^ fi <T M , N ) = (1 - o 1/4 ■ 2 (2) 
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where t = (sinh2E v /k B T sinh2£ ,/l /k B T) 2 , and where f™ N are n-fold integrals. 

The form factor expansions ([1} and (f2]i are of great importance for the study of the 
magnetic susceptibility of the Ising model 

^ = FT ' Yj{(°-o,oo- m ,n) - M 2 }, (3) 

where M = (1 - f) 1 ' 8 f° r T < T c and equals zero for T > T c is the spontaneous 
magnetization. The study of this susceptibility has been the outstanding problem in the field 
for almost 60 years. The susceptibility is expressed in terms of the form factor expansion as 

k B T- X (T) = {l-tf 14 -YjX^T), (4) 



where 



x (n) 



with m = 2n, for T < T c , and m = In + 1, for T > T c , In the last twelve years a large number 
of remarkable properties have been obtained for both x ( "\T) [4|-[ 13 1 and the specialization 
to the diagonal [14| 

Xf(t) - 2 .C- (6) 

These remarkable properties of x {>i) an d Xj(f) must originate in properties of the fj*' N 
themselves. 

For 40 years after the first computations of Wu, the form factor integrals for n > 2 
appeared to be intractable in the sense that they could not be expressed in terms of previously 
known special functions. However, in 2007 this intractability was shown to be false when 
Boukraa et al [ 15 1 discovered by means of differential algebra computations on Maple, using 
the form for the form factors proven in [fl6l , many examples for n as large as nine that the 
form factors in the isotropic case E h - E v can be written as sums of products of the complete 
elliptic integrals K(t 1 ^ 2 ) and E(t 1 ^ 2 ) with polynomial coefficients, where for the diagonal case 
(M = N) we may allow E v + E h . 

These computer derived examples lead to the obvious 

Conjecture 1 

All n-fold form factor integrals for Ising correlations may be expressed in terms of sums 
of products of one dimensional integrals with polynomial coefficients. 

The first discovery that the n-fold multiple integrals which arise in the study of integrable 
models can be decomposed into sums of products of one dimensional integrals (or sums) was 
made for the correlation functions of the XXZ spin chain 

oo 

H XXZ = - J] + «r;«r;., + A,r,r M |. (7) 

7=-oo 

These correlations were expressed as multiple integrals for the massive regime (A < -1) 
in 1992 JT7) and in the massless regime (-1 < A < 1) in 1996 US). In 2001 Boos and 
Korepin |fl9l discovered that for the case A = -1, the special correlation function (called the 
emptiness probability) 



p(n)= <n -r 1 

7=1 V 



>, (8) 



for n — 4 could be expressed in terms of £(3), £(5), £ 2 (3) and In 2, and this decomposition 
in terms of sums of products of zeta functions of odd argument was extended to P(5) in [20] 
and P(6) in ll2~D . Similar decompositions of the correlation function {cr^cr^) were obtained 
for n = 3 in 11221 . for n — 4 in Il23l and for n - 5 in [24|. The extension to the XXZ model 
chain (0 with A + -1 of the decomposition of the integrals for the third neighbor correlation 
(o-' <t' 3 ) for i = x,z was made in ||25l . 

The discovery in |fT31 that a similar reduction takes place for Ising correlations thus leads 
to the more far reaching 

Conjecture 2 

All multiple integral representations of correlations and form factors in all integrable 
models can be reduced to sums of products of one dimensional integrals. 

If correct this conjecture must rest upon a very deep and universal property of integrable 
models. 

In lfT31 the form factors were reduced to sums of products of the complete elliptic 
integrals K(t 1 ^ 2 ) and E(t 1 ^ 2 ). However, the results become much more simple and elegant 
when expressed in terms of the hypergeometric functions F^ and Fjy+i where 

F N = 2 F 1 ([l/2,N+l/2]; [N + I]; t) (9) 

appears in the form factor for n — 1 



JV/2 r \ 



where 

(1/2)a 



(10) 



(11) 



AH 

and (a)o = 1 and for n > 1 (a) n = a(a + 1) • • • (a + n - 1) is Pochhammer's symbol. Note that 
F = fat 1 ' 2 ) = /<>)■ 

The expressions for f^\i{f) in terms of F^ and Fn+\ are obtained from 0151 . rewritten 
by use of the contiguous relations for hypergeometric functions, and we give some of these 
expressions in Appendix A In all cases studied the form factors have the form 

n-\ In 

/<*>) = J] ■ ffft\t) + 2 C™(N; t) ■ F 2 N - m ■ F% +1 , (12) 

m=0 m=0 
A2n+ 1 ) , . v n-1 /<2m+ 1 ) ^ 2n+ 1 

= Z^-^T 1 ^C^\N;t)-F^.F^, (13) 

m=0 ;n=() 

where = 1. The degrees of the polynomials Cm (N; t) are for N > 1 

dBgdj?\N;t) = degCi 2 ' ,+1) (^V;f) = n ■ (2N + I), (14) 

with C, ( ,"'(Af; t) ~ t m as f ~ 0. 

These polynomials are different from the corresponding polynomials in the K, E basis in 
that they have the palindromic property 

C% n \N;t) = f" (2A,+l)+m • C% a \N;l/t), (15) 
C <2n+1) (N;t) = f (2N+l)+m ■ C% n+l) (N;l/t). (16) 
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We conjecture that these results are true generally. 

In this paper we begin the analytic proof of Conjecture 1 and the derivation and 
generalization of the results of |fT31 for the diagonal correlation M — N by studying the 
three lowest order integrals f^ N {t) for n = 2, 3, 4. The results are summarized in Sec. 2. 

In Sec. 3 we derive the results for fj? N (t). We proceed by first differentiating the integral 

f^ N (t) with respect to t, which removes the term proportional to fL jL(Y) from the general form 
( IT2b . The resulting two dimensional integral is then seen to factorize into a sum of products of 
one dimensional integrals. This factorized result is then compared with the derivative of (fT2l 
to give three coupled first order inhomogeneous equations for the three polynomials C m (N; t). 
These equations are decoupled to give inhomogeneous equations of degree three which are 

(2) 

explicitly solved to find the unique polynomial solutions C m (N; t). 

In Sec. 4 we extend this method to fS]f(f). The first step is to apply to f2l,(i) the 
second order operator which annihilates fLlf(i). However, in this case we have not found 
the mechanism which factorizes the resulting three dimensional integral. Instead we use the 
property discovered in ITT31 that the resulting integral satisfies a fourth order homogeneous 
equation which is homomorphic to the symmetric cube of a second order operator and thus a 
factorized form is obtained. This form is then compared with the form obtained by applying 
the second order operator to the form (fi"3) . and from this comparison we obtain 4 coupled 
inhomogeneous equations for the 4 polynomials C^\f). These equations are then decoupled 
to give inhomogeneous equations of degree 5 for C 3 (N; t) and of degree 8 for the three 
remaining polynomials. We then solve these equations under the assumption that a polynomial 
solution exists. 

The results for f^ N {t) with n = 1,2,3 have a great deal of structure which can be 
generalized to arbitrary arbitrary n. Of particular interest is the fact that fff$(f) vanishes 

as f {N+n) and f^" N +1) (t)/t N/2 vanishes as t n{N+n+l) at t -> while each individual term in the 
expansions Sl7\i and ( fT3l vanishes with a power (which may be zero) which is independent 
of N. This cancellation for /Stiff) and fSlf(t) is demonstrated in Sec. 5 and gives an 
interpretation of several features of the results obtained in Sees. 3 and 4. It also provides an 
alternative form (113 81 > for fS N (t) compared to the form (fT3l . In Sec. 6, in a differential algebra 
viewpoint, the canonical link between the 20-th order ODEs associated with the C^{N; t) of 
fxlf(t) and the theory of elliptic curves is made very explicit with the emergence of direct 
sums of differential operators homomorphic to symmetric powers of a second order operator 
associated with elliptic curves, and in an analytical viewpoint, is made very explicit with 
exact expressions (given in Appendix G), for the polynomials C^p(N; t) , valid for any N. We 
conclude in Sec. 7 with a discussion of possible generalizations of our results. 



2. Summary of formalism and results 

The form factor integrals for the diagonal correlations are Ifl5i ri6l for T < T c 



f (2n) (f) = 
Jn,n w 



n(N+n) 



(«!) 2 7T 2 « 



1 2/7 n 
k=l ;=l 



( (l-tx 2j )(x-j-l) \ 



1/2 



(1 - tx 2j -i)(x 2 ] j l - 1)J 

n n ( i-ttLnJ n ^-i-wfe,-^) 2 , 

!</'<« \<k<n v ^ ' \<j<k<n 



(17) 
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and for T > T r 



Jn,n W 



(n+l/2)N+n(n+V) 



nl(n + 1)!tt 2 " +1 



J^l 2n+l n+l 
i ;- 1 



n«i(.-« Sj ) «]-!)]■« n n (r^J 

=1 l<j<n+l \<k<n v 7 ' 

\\ (X2j-X 2k ) 2 - 



n ^- 

1 < j<k<n+ 1 



(18) 



1 < )<£<« 



When f = the integrals in (fTTT i and (fT8l reduce to a special case of the Selberg 
integral ||26l|27l 



f (2n) (f) 
Jn,n"> 



t n(N+n) 



r(Ar + n+l/2)r(/z+l/2) 
(n!) 2 ^ 2 " r(AT+ l/2)r(l/2) 



7=0 



r(jy + ; + i/2)ro- + i/2)rQ- + 2) 

r(;V + n + ; + 1) 



and 



f {2n+1) (t) 



f iV(« + i/2)+«(«+i) r(A f + i/2)r(l/2) 



n!7r 

n-l 

./=« 



! 77-21+1 



r(/V + n + 1) 



r(/v + j + 3/2)r(j + 3/2)ro- + 2) 
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T(N + n + j + 2) 



In particular 



f (2) (?) 

f <3) (f) 



+ i 



i2 

A Af+l 

(2/V+l) 



+ 0(^ +2 ), 



f 3N/2+2 



2(2N + \){N + 2) 2 



+ 0(t 



3/V/2+3x 



(19) 



(20) 

(21) 
(22) 



2.1. General Formalism 

For the special case fSif(t) we will analytically derive the form (fT2l without making any 
assumptions. However, for the general case we will proceed by assuming the forms ( TTZb 
and ( fT3l ) as an ansatz and with this as a conjecture, we will derive inhomogeneous Fuchsian 
equations for the polynomials C„ (N; t) 



Q.^(N-t)-C^(N;t) = tXN;t), 



(23) 



where Q&\N;t) is a linear differential operator and I„ (N; t) a polynomial. 

In all cases which have been studied, the operator £l„(N; t), corresponding to the lhs of 
( f23l , has a direct sum decomposition where each term in the direct sum is homomorphic to 
a either a symmetric power or a symmetric product for different values of N, of the second 
order operator 



, 1 + N — Nt 4 + AN —t- 2Nt 



(24) 
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where D, = d/dt. The operator 2 (N;t) is equivalent to the operator L 2 (N;f) which 
annihilates fLlf(i) iTTBI , as can be seen in the operator isomorphism 

2 (N;t)-t N/2+1 = f /2+1 -Li(N\t). (25) 
The solutions of 2 (N; t) are expressed in terms of hypergeometric functions by noting 

that 

t 2 ■ (1 - t) ■ 2 (N) = t ■ (tD, + a)(tD, + b) - (tD, - a')(tD, - b'), 

with 

a = -A/ -1/2, b = -1/2, a' = N + l, b' = 1, (26) 

which for \t\ < 10 has the two fundamental solutions Il28l p. 283] 

f' ■ 2 F l ([a + a',b + a'];[a'-b' + l];t), t y ■ 2 F l ([a + b',b + b']; [b' - a' + 1]; f). (27) 

Using d26l i we have the two solutions of 2 (N) 

Ul (N;t) = f A ' +1 • 2 Fi([l/2, 1/2 + AH; [N + 1]; f) = t N+l ■ F N , (28) 
and: t ■ 2 F 1 ([l/2, 1/2 - N]; [1 - N]; t). (29) 

The solution U\(N\ t) in d28b is regular at t = and has the expansion 

CO 

u i (N; f) — f w+I • Y b n (N) ■ t", (30) 

with 

„„ m . W ?f> 2 *^ . (3D 
(N+ 1)„«! 

Since we will in this paper work with positive integer values of N, it is better to introduce 
as the second solution 

t N+1 ■ 2 ^i([l/2, 1/2 + AH; [1]; 1 - f). (32) 



When Af is not an integer the hypergeometric function (1321 can be written as the following 
linear combination of the two previous solutions (l28t and (l29l 

If ~ A) f +l ■ 2 Fi([l/2,l/2 + AT|; [N+l]; t) 



r(i/2)r(i/2-A0 
r(A0 



r(l/2)T(l/2 + A0 



■t ■ 2 ^i([l/2, 1/2 - N]; [l-N]\ t). (33) 

The hypergeometric function (132) is not analytic at t = but, instead, has a logarithmic 
singularity. 

X For |t| > 1, we write z= 1/f and the identical procedure is found to interchange a with a' and with //. Thus the 
two fundamental solutions valid near t = oo are Ui(N;z) = z~ l/2 • %F i ([1/2, 1/2 + AT]; [1 + = ?r >/2 ■ F N , 

U2(N; z) = zr N ~^ 2 ■ iF\ ([1/2, 1/2 — JV]; [1 — iV]; z). The identification of the hypergeometric functions of <28) and 
{29} with these two solutions is a consequence of the palindromic property of the operator Oi(N; t). However, we 
note that Uj(N; z) is not the analytic continuation of uj(N\ t). 



From Il29l (2) on p. 74 and (7) on p. 75] we may choose to normalize the analytical part 
of the second solution to t as t — > 0. Denoting such a solution U2(N; t), it reads 



JV-l co 

u 2 {N- 1) = t ■ J] a„(N) ■ f + t N+x ■ N ■ 4 • J] b n(N) [k„ - ln(0] ■ f, (34) 

n=0 n=0 

with a (N) = 1 and for n > 1 

(l/2)„(l/2-jV)„ (l/2)„(jV-n)! 
(1-A0„«! (l/2)iv_„n! 

and k„ = H B (l) + H„ +N (l)-H„(l/2)-H n+N (l/2), where 

n-l 

z ■ 



n-l j 

h z+k 



(36) 



are the partial sums of the harmonic series. The series expansion ( f34T > corresponds to the 
maximal unipotent monodromy structure of C>2(N; t) which amounts to writing the second 
solution as: 

u 2 (N;t) = w 2 (N-t) - N-A 2 N - Ul (N-t)-\n(t) (37) 

where w 2 (N; f) — t + ■ • • is analytical at t - 0. This function w 2 (N; t) is the solution 
analytic at t =0, different from uy(N; t), of an order-four operator which factorizes as the 
product 2 (N; t) ■ 2 (N; t), where 2 (N; t) and 2 (N; t) are two order-two homomorphic 
operators 

d 2 (N;t)-h = h-0 2 {N;t), (38) 
where one of the two order-one intertwinners I\ and J\ is quite simple, namely 

1 t -2 N 

h - r D -~2-fir^)-2f- (39) 

Finally, we note the relation which follows from the Wronskian of 2 (N; t), 

Mi (AO • u 2 (N + 1) - p N ■ u 2 (N) ■ ui(N + 1) = f +2 , (40) 

with 

(2N+ l) 2 

B N = — . (41) 

H 4N(N +1) 

( 2 ) 

2.2. Explicit results for f NN {f) 

For f^ N (t) the parameter and the polynomials (N; t) of the form (fT2l i are explicitly 
computed in Section[3]as 



and 



K {2) = N/2, (42) 



2;V+ 1 -m 



C%\N;t) = A^-r- J] c%{N)-f, (43) 



B=0 
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with 

Af = (-ir 1 -^-^.^. (44) 

Using the notation that 

lf]n = = the coefficient of f in the expansion of / at t — (45) 
we have for < n < N - 1 

n 

eg(A0 = c ( ® 2N _ 1 _ n {W= [r 2 u 2 (N) 2 ] n = J] flt (A0 • fl„-*(A0, (46) 

n 

= C i?w-»( JV )= [r 2 M2(^V)«2(^V+l)]„ = 2^(A0-a„-*(iV+l), (47) 

yt=0 

and 

'l;AA J V ~ " TL 2,A<-1 

and where for < n < N 



cfhN) = 4 +c® AN), (48) 



c o?„W = = [r 2 u 2 2 (N+l)] n = c^ B (iV+l), (49) 

where a n (N) is given by (f35T >. We note that the sum (l46l ) for eljL _j may be written by use of 
the second form of a n (N) in ( f35l > in the alternative form 

eg,., = 4- 2N-H N (l/2), (50) 

where ffjv(z) is given by ( l36l >. 

We also derive the recursion relation for N > 1 

0)= -i^-z ^ ^TiT (51) 

2.3. Explicit results for f N N (t) 

For fSlf(t) the parameter Xq 3) and the polynomials C^iN; t) of the form (IT3l are explicitly 
computed in Section|4]as 



K f = — — - (52) 
and 

2N+\-m . 

C^iV; = A® • r • Yu c >n,n(N) ■ f + -j=-^ • Cg>(JV, 0, (53) 
where we make the definition C?\N, t) — and 



3N+1 



= (-lr 1 -!-^)-*^. (54) 



The coefficients c^ n (N)'s are given by a simple quartic expression of the a„'s and b n 's. For 
< n < N - 1 they read 

= 4;W- 2 -„W = [r A, - 4 ^(^v)« I (^v)]„ = 

n m I 

= J] Z Z ak(N) ' a '- k(N) ' fl '"- /(A0 ' ^-"^ ^ 55 > 



and 



m=0 /=() £=0 



« m I 

= Z Z Z flt(A ° ' fl " (A ° ' fl, "-' (A? + 1} ' (56) 



m=0 /=() fc=0 

for < n < N 



m=0 /=() <:=() 

and for < n < N - 1 



Z Z Z fli(A? + 1} ■ a/ -* (A? + 1} ■ flm -' (A? + 1} • b "-"^ (5 ^ 



n m I 



= Z Z Z °* (A ° ■ a " (A? + 1} ■ fl *- m(A? + 1} ■ < 58 > 



=0 /=0 jfc=0 



with the middle term of C? \N; t) of order N + 1 



c%= ^ + ^-^-[(iV-l). e&_ 1+ 4e&_J 

|l^'K-4^^-l/4),y, (59) 
where a„(2V) and b„(N) are given by (l35l > and (fJTJ. 

3. The derivation of the results for f ( ^ } N (t) 

We begin our derivation of the results for fj®, of Sec. I2.2l bv integrating (fTTT i (with 2n = 2) 
by parts using 

u = / r - 1/2 .(l-y) 1/2 .(l-0') 1/2 , (60) 
, [JV ■ (1 - y) ■ (1 - fy) - 1/2 (1 -fy 2 )] 

rfM = y (l-v)'/ 2 -(l- f v)'/ 2 dy > m 

dv = (Tw' v = r^' (62) 
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to find 



f (2> (t) 



rl pi ( N+l 

Jo ^ Jo ^ 2? 



^+1/2^-1/2 (1-^) 



2tt 2 (1 - j) 1 ' 2 (1 -^^(l-y) 1 / 2 ^ - fy) 1 / 2 (1 - txy) 

ri f N + l x N + M2 y N-ll2 {l _ y) ll2 {l _ ty) ll2 



r r t 

■N dx dy—r- 
Jo Jo 



, jo n- (l-x) 1 /2(l- fx) l/2(l_ f ^) 

The first term in (l63l is separated into two parts as 

f 1 f 1 f-x^x 1 / 2 1 

m ay 
Jo Jo 



(63) 



27T 2 y/ 2 (1 - x yi 2 (1 - fx) 1 / 2 (1 - y) 1 / 2 (1 - fy) 1 / 2 (1 - fxy) 



^ 2tt 2 (l-x) 1 / 2 (l-fx) 1 / 2 (l-y) 1 / 2 (l-fy) 1 / 2 (l-fxy)' 1 ^ 



- \ dx \ 

Jo Jo 



and in this second term we interchange x <-> y. Then, recombining these two terms, we see 
that the factor 1 - txy cancels between the numerator and denominator in (164V Thus the first 
term in (|63l factorizes and we find 

JM pi WV+1 ^+1/2^-1/2 
</v </v 
o Jo 



2n 2 (1 -x)'/ 2 (l -fx) [ / 2 (l -y)'/ 2 (l - fy) 1 / 2 

' ri ,Af+l x ^+l/2 3; iV-l/2 (1 _ y) l/2 (1 _ f3 ;)l/2 



ofx I dy — r- 

Jo K 2 



(1 -X) [ / 2 (1 -fx)!/ 2 (l -fXj) 

= I f (l) f (l) 

2 ->N,N Jn+i,n+i 

C 1 C 1 t N+l x N+ll2 y N -^ 2 (\-yyi 2 {\-tyy^ 2 

- N Jo dx Jo dy — T?W^F' (65) 

From 07J we find for N > 1 that the integral in the second term of d65j is - /®j JV+1 (f) 

and thus we have 

A/2 

O') = — •/^(o-4ViW -^-O f ) -C, w+I (0]- (66) 

From (|66T > we obtain the recursion relation 



A? + 1 „™ . . f 1/2 

-W') = 

and thus for N > 1 



A2) ( ) fvj^i . ,(2) ( } _ f_ < ( } _ ( ) (6?) 



Oo= -f^ 1/2 -Z - (j T!) ■ (68) 



To proceed further we return to (163) which we write in terms of F N as 

,(2) , A _ jjV ^jV+1 JV+l p p 



- N f dx f 
Jo Jo 



1 f N + l x N+l/2 y N-l/2 (l _ y y/2 il _ ty A/2 

dy ~W (1 - x) 1 ' 2 (1 - fx) 1 / 2 (1 - txy) ' (69) 
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The integral in (l69l does not have a manifest factorization. However, if we compute 
df!?L(f)/dt in the contour integral form of (IPTl , and note that 



d 
dt 



(y-f 1/2 )(i-f 1/2 y) 



(x-t l ' 2 )(l -t l > 2 x) 



1/2 



1 

Jjl 



(y 



(x-r'/ 2 )(l -t^ 2 x) 



1/2 



(jcy-l)U-y)(f-l) 



(70) 



(y - f!/ 2 )(l - t^ 2 y)(x - t^ 2 )(l - t l ' 2 x) ' 

the resulting integral does factorize and, introducing Gn, some well-suited linear combination 
of F N and F N+1 , 

G N = 2 F l ([3/2,N + 3/2]; [iV+1]; f) 

1 +f f 2N + 1 



we find 

dt 



(1 



(1-f) 2 



(2jV + l)4 
16 (A? +1) 



(1-f) 2 N+l 



Af+1> 



(71) 



[(2iV+ l) 2 -^ + i -Gjv 

~(2N-l)(2N + 3)-F N -G N+l ]. 



(72) 



It remains to integrate ((72). However, in general, integrals of products of two 
hypergeometric functions with respect to the argument will not have the form of the product 
of two hypergeometric functions. We will thus proceed in the opposite direction by 
differentiating (fT2l for 2n — 2 with respect to t and equating the result to (l72l to obtain 

(2) 

differential equations for the C m (N; t) which we will then solve to obtain the final results 

From a straightforward use of the contiguous relations of hypergeometric functions [29 1, 
we introduce the following well-suited linear combination of F^ and F^ + i 

4(N+1) dF N 

F N = 2 F l ([3/2,N + 3/2]; [N + 2]; t) = 



2N+1 dt 
= _L. (2.(^+1).^ -(2N+1)-F N+1 ). 
The derivative of (fT2l with 2n — 2 may be written in the quadratic forrrS 



B\ ■ + Bi ■ F N ■ F N + B3 • F 2 ,, 



(73) 



(74) 



with 



dC, 



(2) 



dt 



(N+l) 
2(W+l/2)f 



■Cf + 



(N+l) dC 



(2) 



(N+l 12) dt 



(N+l) 2 
(N+l/2) 2 t 



c 



(2) 



(N+\) 2 dC 



(2) 



(N+l/2) 2 dt 



(75) 



B 2 = 



(N+l/2) 
(N+l) 

1 + 



r (2) 
u 



1 



1 -2/ + AT(l -0 



2(N+l/2) 

(N+l)(3+2N-2t) 
2(JV+l/2) 2 f 



2(N+l/2)t 



C 



(2) 



(1 - f) dC 



C 



(2) 



(N+l)(l-t) dC\ 



(2) 



(JV+1/2) 2 dt 



(2) 



2(JV+l/2) dt 



(76) 



§ For convenience the dependence of the C$ on N and / is suppressed here and below (see J781 ). 
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n - (1 ~ f) r® - gjj^lfXlzg r P) . d- f ) 2 flfC 2 2) r77 , 

3 4(7V+1) 2 4(JV+l/2) 2 r 2 4(JV+l/2) 2 A ' K ' 

The derivative of fSit(t) in d72l by use of contiguous relations |29| is expressed in terms of 
Fn and as 

<O0 



where, 



7?4 - -. - AjU^j 



N,N^ ' 2 ~ ~2 

= Z? 4 ■ Fjy + B5 ■ F N ■ F N + Z?6 ' (^S) 

2N + 1 

4 

[f -jV(l-Q](2jV+l) j2 „ 

B5 ~ WTT) v ' 

6 4(JV+1) 

(2) 

5.7. Linear differential equations for C m (N; t) 

To obtain the C ( ^(N; t) we equate d74l i with (TTBl and find the following first order system of 

(2) 

equations for C m (N; t) 

(2N+1) 2 „ £Cg (/V+l) (2) (jV+1) rfCf 
4 ' N ' dt 2(N+l/2)t' 1 (N+l/2) ' dt 

(N+D 2 M2) , (/V+l) 2 «*Cf 



Q + ,„ : , ■ — r-, (80) 



(N+l/2) 2 ? 1 (N+l/2) 1 A 



(2N + l)-[t-N(l~t)] , 2 „ (iV+1/2) 



,2 JV UV + V^ ~(2) 

/Ixf * ( — ' La 

4(tf+l) w (AT+1) 

-(2) 

-, , [ I — / ) < f v_ 



j + i ^ i-2r +jV(i-o l ^ (2) (i - o rfq 4 - 



2(7V+l/2) 2(7V+l/2)f 



' l 



2 (N+l/2) dt 



(N+l)(3+2N-2t) (2) (JV+Dd-Q dCf 
2(7V+l/2) 2 f ' 2 (7V+1/2) 2 dt ' 



(2iV+1)2 .4.(1-0.^ = -^L-cf> 



16(7V+1) 2 ™ 4(7V+1) 

(2 + 2JV-Q(l-0 | (1-Q 2 </Cf 

4(7V+l/2) 2 f 2 4 (N+l/2) 2 dt 

From this first order coupled system we obtain third order uncoupled equations for the 

Ci 2) (/V;r) 

2(l-f) 2 -f 2 -^f- -6(N - (N -l)r)(l -t)t--^- 

dC {2) 

+ 2[N + 2N 2 + (1 + 4iV - 47V 2 ) t -(5N- 2N 2 ) t 2 ] — 

dt 

m N(N+1)(2N+1) 2 , v 
+ (2N + l)(27Vf-27V- 1) • Q 2) = — ^ --^ •(!-*) -A (83) 
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d 3 C {2) d 2 C {2) 
2(1- tf (1 + t) ■ f 3 • — i- - 2(1 - [1 + 3N + At + (1 - 3N)t 2 ] ■ t 2 ■ — -f- 

<iC (2) 

+ 2 [2 + AN + 2N 2 + (3 + AN - 2N 2 ) • r - (3 + 8tf + 2N 2 ) ■ t 2 + 2 N 2 ■ f 3 ] ■ t — 

dt 

-[A + 8N + AN 2 + (5 + 6tf) • r - (5 + ION + AN 2 ) ■ t 2 ] ■ cf ) 
\2 . r i \t2 t\r _i_ . ft _i_ ^ 



(2jV + l) z ■ [-2/^ (jV + 1) ■ (r + l) 2 + (4tf + 1) ■ t] 2 „ +I 

and 



2(1- f) 2 • f 3 • — rf- - 6 (1 + - M) • (1 - t) ■ t 1 ■ — -f- 

dr atr 

dC (2> 

+ 2 [7 + 9N + 2N 2 - (7 + 127Y + 4tf 2 ) • t + (1 + 3tf + 2tf 2 ) • f 2 ] • t ■ 2 



off 

[16 + 2AN + 8N 2 - (15 + 28tf + 12N 2 ) • t + (2 + 6N + AN 2 ) ■ t 2 ] ■ C (2) 
N 2 (2N+1) 4 - (1 -t) 



8(N + l) 2 



•4-^ +2 - (85) 



From d83l-<f85t it follows that C (2) (N;t) and t 2N+m+l • C (2) (N;l/t) satisfy the same 

(2) 

equation and thus, if C m (N;t) are polynomials they will satisfy the palindromic property 
JBV From dS3) and (H) it follows that the polynomials C (2) (N ; f) and C 2 2) (tf ; f) satisfy 



C, (2) (tf;f) = = ■ C (2 \N +\;t). (86) 

(N+l)-(P N+1 -fi 2V 



We therefore may restrict our considerations to Cf\N; t) and C^\N; t). 

We will obtain the polynomial solutions for the differential equations (l83l-(l85Tl by 
demonstrating that the homogeneous parts of the equations are homomorphic to symmetric 
products or symmetric powers of the second order operator 02(N). 

3.2. Polynomial solution for C? (N; t) 

Denote ctf?(N,i) the order-three linear differential operator acting on cf\N,t) on the left 
hand side of dSBI l. Then it is easy to discover that the operator Q^(N, t) is exactly the 
symmetric square of the second-order operator OziN; f) 

Q. i2 \Nt) = Sym 2 (0 2 (tf ; f)), (87) 
which has the three linearly independent solutions 

u x {N;t) 2 , ui(N;t)-u 2 (N;t), u\(N;t) (88) 

where the functions uj(N; t) for j = 1, 2 are defined by (l30ll-(l3"6il. The indicial exponents of 
d85l ) at t = are 

2N + 2, N + 2, 2, (89) 

which are the exponents respectively of the three solutions fl88l >. Therefore, because the 
inhomogeneous term in d85l l starts at t N+1 the coefficients c^ n in d43l for < n < N — 1 
will be proportional to the first N coefficients in the expansion of u\(N ; f) about t — 0. 
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Equation (|85l l is invariant under the substitution 

Cf{N;t) — > f 2A,+3 -Cf (A^;l/f), (90) 

which maps one solution into another. Therefore if it is known that the solution C 2 (N; t) is a 
polynomial the palindromic property 

JV - J® (91 ) 

must hold and thus C 2 \N; t) is given by d46b where the normalizing constant A!, 2 ' remains to 
be determined. 

However, the invariance ( |90l is by itself is not sufficient to guarantee the existence of 
a polynomial solution with the palindromic property (PT3V To demonstrate that there is a 
polynomial solution we examine the recursion relation which follows from d85l l 

Af -{2n(2N- ») (N- n)-cg (N) 

+ (ANn - 2N - 2n 2 + 2n - 1) (2n - 1 - 2A0)) • c®_,(iV) 
+ 2 (n - 1) (2N- n + l)(N - n + 1)) • cg_ 2 (AT)} 
jV 2 (2jV+l) 4 2 

= (0„,N - 0„Jf+l) ■ g(jv+ • Ajy. (92) 

where c 2 . n (N) = for n < -1 and we may set c!,.q = 1 by convention . By sending 

n — » 2A^ - n + 1 in ( |92] l we see that c^ n {N) and c^iv-n-i^) ^° sat i s fy me same equation as 
required by (|9TV 

(2) 

To prove that the solution C 2 (N; t) is indeed a polynomial we examine the recursion 
relation d9"2i l for n — N. If there were no inhomogeneous term then, because of the factor 
jV - n in front of c 2 2 ) , the recursion relation (l92t for n — N would give a constraint on c ( 2 Z> N _ l 

and c 2 2) N _ 2 . This constraint does in fact not hold, which is the reason that the solution u 2 (N; t) 

JV+2i 



is not analytic at t — but instead has a term t In t. However, when there is a nonzero 
inhoi 
Thei 



inhomogeneous term at order t N+1 the recursion equation (l92i l is satisfied with a nonzero A 2 . 

(2) 

The remaining coefficients c 2 ' for < 2N - 1 are determined by the palindromy constraint 



(2) (2) 

For C 2 (N; t) to be a polynomial we must have c\.'(N) = for n > 2N. From the 
recursion relation (|92l we see that because of the coefficient 2N - n in front of c^(A^) the 

(2) (2) 

coefficient c 2 . 2N (N) may be freely chosen. The choice of c 2 . 2N (N) + corresponds to the 

(2) 

solution of £l\ (N; t) which has the indicial exponent N + 2 and clearly does not give a 
polynomial solution. However by setting n = 2N + 1 in ( |92] l we obtain 

2{N+\){2N+\)-c% N+l (N) - (2N+\) 2 -c% N {N) = 0. (93) 

and if we choose c^L^AO = Owe obtain c ( 2 l 2N+ - l (N) = also. Therefore because (l92l is a three 

term relation, it follows that c?}(N) = for n > 2N as required for a polynomial solution. 

It remains to explicitly evaluate the normalization constant A*? which satisfies d92| > with 
n — N. A more efficient derivation is obtained if we return to the original inhomogeneous 
equation d85l l. Then we note that if we include the term with n = in the second terms on 
the right-hand side of (l34l l in the computation of the term of order f A ' +2 in the left hand side 
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of d85l l we must get zero because w 2 is a solution of the homogeneous part of ( |85l l. Therefore 
when we use the extra term in u\ of 

- 2 f +2 ■ N ■ A 2 N • In f, (94) 
in the lhs of (l85l . and keep the terms which do not involve In f, we find 

2(N 2 -l)-c ( 2 %_ 2 (N) -(2N 2 -l)-c™ jtfV) = -4^V 3 -4. (95) 



Thus, using d95T > we evaluate d92l ) with n = N as 

1(2)^3 3 2 



jV 2 (2JV+l) 4 j2 
8(JV+ l) 2 



-4A^W J 4 = \ ' (96) 



and thus 



Af = (97) 



3.5. Polynomial solution for C\ (N; t) 

The computation of cf\N; t) has features which are characteristic of C%\N; t) which are not 

(2) 

seen in C 2 (N; t). Similarly to what has been done in the previous subsection we introduce 
Q.f\N; t), the order-three linear differential operator acting on cf\N; t) in the lhs of (|84| >. 

(2) 

The indicial exponents at t — of the operator Qj (N; t) are 



1, N+l, 2N + 2 (98) 
) is found to b 

and 02(N + 1) by the direct sum decomposition 



(2) 

This order-three operator Qj '(N; t) is found to be related to the symmetric product of C>2(N) 



Sym(o 2 (A0, 2 (N + 1)) • t = Q (2) ©(d, (99) 

The three linearly independent solutions of £l^(N; t) are to be found in the set of four 
functions 

f 1 ■ u\{N; t) ■ ui(N + l;f), f _1 • u 2 (N,t) ■ ui(N + l;t), 

r 1 ■ ui(N;t) ■ u 2 (N+ l;f), f" 1 ■ u 2 (N; t) ■ u 2 (N + 1; t), (100) 

where from the definitions of U\(N; f) in ( f30b and M2(^V; f) in d34l > the behaviors of these four 
solutions as t — > are f 2iv+2 , f^ 2 , f^" 1 " , f respectively. 

Following the argument given above for C 2 (N; t) we conclude that because the 
inhomogeneous term in ([83} is of order t N+l that the terms up through order must be 
proportional to the solution of the homogeneous equation 

r l ■ u 2 (N;t) ■ u 2 (N+ 1;0, (101) 

which begins at order t. This observation determines the form ( 143) and the coefficients d47| > 
cf}(N) for < n < N - 1 . The normalizing constant A (2) and the remaining coefficient cf^ N (N) 

(l48b are then obtained from the inhomogeneous equation (1831 . Finally, to prove that cf\N; t) 
is actually a palindromic polynomial the recursion relation for the coefficients (AO must be 



used. Details of these computations are given in Appendix B 
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3.4. The constant K {) 



(2) 



Finally, we need to evaluate the constant of integration K Q in (fT2l >. This is easily done by 
noting that from the original integral expression ([PTt that f^Li®) = for all N. From (l43Tl- 
( 144b we see that 

C|, 2) (^;0) = --, Cf\N;Q) = cf{N;Q) = 0, (102) 

and using this in ( fT2l we obtain K ( { 2) = N/2 as desired. 

4. The derivation of the results for f$l/(t) 

The form factor fSli(f) is defined by the integral (TT81 with 2n+l - 3, and if we are to follow 
the method of evaluation developed for f^ N (t), we need to demonstrate analytically that there 
is an operator which, when acting on the integral, will split it into three factors. Unfortunately 
we have not analytically obtained such a result. 

However, we are able to proceed by using the methods of differential algebra and 
from [15 1 it is known computationally for integer N that fjfy is annihilated by the operator 
U(N) ■ L 2 (N) where 

, 2t - 1 1 1 N 2 

L 2 (N) = D 2 + D, + -, (103) 

y ' ' (t-l)t At 4(f-l) 4f 2 v ' 

and Lq,(N) annihilates f ( ^ N (t), and where, 

4 (2f - 1) , (241f 2 -241f + 46) , 



(2f- l)(122f 2 - 122f + 9) 81 (5f - l)(5?-4) 5 N 2 

(/-l) 3 / 3 °' + 16 f 3 (f-l) 3 2~ 

(23 -32t)N 2 9 (8 - 17 t)N 2 9N 4 



+ - -^-D, + -- —— + -. (104) 

2(/-l)f 3 8 (f-l)f 4 16 f 4 

Furthermore the operator Ln(N) is homomorphic to the symmetric cube of L2(N) by the 
following relation, 

L 4 (N)-Q(N) = R(N) ■ Sym 3 (L 2 (A0), (105) 

where, 



Q(N) = (r-l).r.D? + l(2t-l)-D 2 + ( ~4lt + 6) 
^ ' J '2 ' 4(f-l)f 



9(2f-l) 9(t-l)N 2 9(2f-l) , 

+ 8(Tr/t-4-^- D '-s--l 2 - 2N > (106) 



and 



, 23 , 21 6 -29? + 29f 2 

R(N) = (t-D-t-B} + —(2t-l)-D 2 +— — — A 

9 (2f-l)(125f 2 -125f+16) 9JV 2 /(f-1) (10f-9) \ 

+ s jj^Wt 2 a ? 2^ (1U7) 
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We therefore conclude that since f2lr(i) i s regular at t — and the solution of L2(N) 
which is regular at t = is F N , that 

Q(N) ■ B ■ t 3N/2 ■ Fl = L 2 {N)-f^ N , (108) 

where Bq is a normalizing constant which is determined from the behavior at t — 0. From 
the integral ( fT8l we find 



(3) 7V + 2 / \ . 3W/2+2 0f< 3ff/Z+3x QO m 



and from the expansion of F^v we have 

8(tV + 1) 2 (tV + 2) 



e(^v)-f 3,v/2 -^ = 3(2jv+D 3 +0(rW/2+ i x (110) 



and thus 



1 



Bo = --4- (HI) 

Operating 2(A0 on t 3N ^ 2 F^, one can write the result in the basis F^ and F^. Similarly, 
one can operate on the form f^L in (fl~3T > with L 2 {N) and write the result in the same basis F N 
and F^. Then, matching powers of the hypergeometric functions on both sides of the relation 
( 1108t will yield four coupled inhomogeneous ODEs to be solved. The four coupled ODEs are 



given in Appendix C 



For Cm'(N; t) with m = 0,1,2, the reduction of the four coupled second order equations 

(3) 

leads to inhomogeneous 8-th order uncoupled ODEs for each C m (N; t) separately, of the form 
Yj P m,j{t)-t j -—.C^\N-,f) = I m (t), (H2) 



=0 



where 



7=0 7=0 

h = l N+2 -Y J hii)-t i , (113) 

7=0 

where the I m (t) are antipalinromic and P m> „(?) are polynomials. In particular 

P m , % {t) = (l-tf-P m {t), (114) 

where Po(t) and P2(t) are order six and P\(i) is order eight. 

However, for C?\N; t) a step-by-step elimination process in the coupled system 
terminates in a fifth order equation instead. We derive and present this 5th order equation 
in '. 



Appendix D but the eighth order equations given by Maple are too long to present. 
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4.1. Polynomial solution for C 3 (N; t) 

IS, of the ODE dDTO for C 3 
to be isomorphic to Sym 4 ((92(A0) • t {N+l \ the symmetric fourth power of OiiN) multiplied 
by f N+l \ Therefore all five solutions of the homogeneous equation are given as r (N+l) times 
products of the solutions u\{N; t) and U2(N; t). The fifth order ODE has at t = the indicial 
exponents 

-N + 3, 3, N + 3, 2N + 3, 3N + 3. (115) 

Therefore because the polynomial solution must by definition be regular at t — the first 
N + I terms (from t 3 through f A ' +3 ) in the solution 

r (N+l) ■ u\(N) -KiCAO, (116) 

which vanishes as f 3 , will solve the inhomogeneous equation ( ID. lb . so that 

2N-2 

Cf{N; t) = Af ■ t 3 ■ Yj 4-» ■ f "' < ! 1? ) 

«=0 

where for < n < iV - 1 

ra m / 

4;« = Z Z Z flt(A?) ' a '- k{N) ' a "'-' (N) ' b "-'" ( W- < 118 ) 

m=0 /=() /t=0 

The lowest order inhomogeneous term is f' v+3 which is the next indicial exponent in ( II 15b and 
therefore the normalizing constant A\ is found from the first logarithmic term in the solution 
of the homogeneous equation by exactly the same argument used for C\ (N; t). Thus we find 



The homogeneous operator on the LHS of the ODE ( ID. It for C?{N, t) is found on Maple 



Af = --/3l-A N . (119) 



2 
3 

The remaining demonstration that C9\N; t) is a palindromic polynomial follows from 
the recursion relation for the coefficients, as was done for Cf(N; t), with the exception that 
because the inhomogeneous term in (ID. It is proportional to t N+l (t 2 - 1) instead of t N+1 (t - 1), 



there is an identity which must be verified. Details are given in Appendix D 



(3) (3) 

4.2. Polynomial solutions for C 2 (N; t) and C\ (N; t). 



(3) (3) 

A new feature appears in the computation of C 2 (N; t) and C\ (N; f). 

(3) 

The indicial exponents at t = of the 8-th order operator £l 2 (N; t) 

-N + 2, 2, 3, JV + 2, JV + 3, 2N + 2, 2N + 3, 3N + 3, (120) 

and for Q. (3) (N; t) are 

-N, 0, 1, N+l, N + 2, 2N + 2, 2N + 3, 3N + 3, (121) 

and from these exponents it might be expected that the solution of Q. ( 2 3) (N; t) (Q.f(N; t)) 
which is of order t 2 (f°) could have a logarithmic term f 3 In t (t In t) which would preclude the 
existence of a polynomial solution of the corresponding inhomogeneous equation. However, 
this does, in fact, not happen because there is a decomposition of the 8-th order operators into 
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a direct sum of the third order operators Q®(iV; t) (Q.f\N; t)) with exponents 2, N + 2, 2N + 
2(0,^+1,2^ + 2) and new fifth order operators M%\N; t) 

Q%\N;t) = M%\N;t) © n%\N;t) (122) 

with exponents -JV+2, 2, AT+2, 2N+2, 3N+2 for M®\N; t) and -A 7 , 0, N+\, 2N+2, 3N+3) 
for Mq \N; t). Furthermore M^\N; t) is homomorphic to the symmetric fourth power 
of 2 (N) and M^\N;t) is homomorphic to the symmetric fourth power of 2 (N + 1) 



(see Appendix E for details). The inhomogeneous equation is solved in terms of a linear 
combination of the solutions of the third order and fifth order homogeneous equations. 
However, a simpler form of the answer results if we notice the isomorphisms 

af(N;t) = Sym(0 2 (N),0 2 (N),0 2 (N),0 2 (N+l))-t N+2 , (123) 
a^(N;t) = Sym(0 2 (N),0 2 (N +l),0 2 (N +l),0 2 (N +l))-t N+ \ (124) 

(3) 

The desired solutions for £l 2 (N; t) are constructed from the two solutions which have 
the exponents 2 and 3, 

r N - 2 ■ uj(N) ■ Ul (N) ■ u 2 (N + i), r N ~ 2 ■ u 3 2 (N) ■ m(N + 1), (125) 

which by use of the Wronskian condition (1401 may be rewritten as a linear combination of 
two solutions each with the exponent of 2 as 

Af ■ r N - 2 ■ ul(N) ■ Ul (N) ■ u 2 (N + I) +Bf-u\(N), (126) 

and similarly for C\'(N; t), we choose as the solution of the homogeneous equation the two 
solutions with exponent 

A® ■ t~ N ~ A ■ u\(N + 1) • wi(A0 +Bft- 2 -u\{N+Y). (127) 

~( 3 ) fr,i- fl <r" *, <r" AT _ 1 QnH onrl 



This procedure determines the constants c\. for < n < N - 1 and and c\. for < n < N. 

c (3) 



with palindromy determining the remaining ci. for N < n < 2N - 1 (l56l and and Cq 3 ' ( for 



AT + 1 < n < 2N + 1 (ED- 

The constants AS, 3 ' and fij 3 ' in ( 1531 are found by using (l53l with d56l in the 

(3) 

inhomogeneous equation for C 2 (N; t) and matching the first two terms in the inhomogeneous 
terms of orders t N+2 and t N+3 (which are the same orders as the corresponding indicial 
exponents (11201 i). This generalizes the determination of A„ for C„ above. Similarly the 
constants Aj 3) and are found using ( l53l with d5"7l in the inhomogeneous equation for Cy 3) 
and matching to the inhomogeneous terms t" and t N+2 . Thus we obtain the results d5"3l-(l54l 
summarized in section |2~3"l 

(3) 

4.3. Polynomial solution for Cj (A 7 ; f) 

(3) 

The computation of Cj (A 7 ; f) has further new features. 

(3) 

The 8-th order homogeneous operator Q t (A 7 ; f) of the inhomogeneous equation for 
Cj (A 7 ; has the eight indicial exponents at t — 



-AT+1, 1, 2, Af+1, N + 2, 2N + 2, 2N + 3, 3N + 3, 



(128) 
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and, as in the case of C)!P(N; t) and (A 7 ; t) has a decomposition into a direct sum of 

(2) 

Qj (JV; f) and a fifth order operator. However, simpler results are obtained by observing that 

f3) 

Qj (N; t) is homomorphic to the symmetric product 

Sym(<9 2 (A0, 2 (N), 2 (N + 1), 2 (N + 1)) • f /v+3 

= Q fV; f) e(A-^f^), (129) 

which satisfies a 9-th order ODE with indicial exponents at t — 

-JV+1, 1, 2, AT + 1, JV + 2, JV + 3, 2JV + 2, 2N + 3, 3JV + 3. (130) 

The solutions with exponents of 1 and 2 are respectively 

f _JV " 3 ■ u 2 {N) ■ u^N) ■ u\{N + 1), r^ -3 -ul(N)-u 2 (N + l)-Mi(iV+ 1), (131) 

and again, recalling the Wronskian relation (|40l , we may construct the polynomial Cf\N; f), 
similar to the construction of cf (N; t), from the linear combination 

Af ■ r N ~ 3 ■ u 2 (N) ■ hi (AO • u\{N + 1) + Bf ■ f 1 • u 2 (N) ■ u 2 (N + 1), (132) 

which determines the coefficients c\. for < n < N — 1, with palindromy determining the 

remaining cf n for N + 1 < n < 2N (gS). The coefficients Af and fif <fl32l are determined in 

a manner similar to the determination of Af and fif , by matching to the terms of order t N+l 
and^ 2 . 

Finally the term c\. N is computed by using the previously determined results for 
Cf (N; and cf \N; t) in the coupled differential equation ( lC.3t , giving the result (|59l 

(3) 

4.4. Determination of 

It remains to determine the constant d52l , which is easily done by setting t = in ( fT3l ) to 
obtain 

= ■ A w + A%> + ^± ■ Aw • Af, (133) 

and using (144-b and (|54l . 

5. The Wronskian cancellation for (f) and f^ N (t) 

The polynomials C„ (N; f) are of order f m as t — > 0. However, from (ETb and (l22l we see 
that fSir(i) vanishes as t N+l and fS N (f) vanishes as t N+1 . Therefore for t — > 0, a great deal 
of cancellation must occur in (TT2l and ( fT3l . This cancellation is an important feature of the 
structure of the results of sec. 2.2 and 2.3. 

To prove the cancellations we note that the «-th power of the Wronskian relation (l40t is 

r" (N+2) ■ is-W ■ (") ■ ^ N + l >iWr j ■ [u 2 (N) Ui (N + l)Y = 1, (134) 

or alternatively, 

7=0 VI L 



1 ■ u 2 {NV ■ F n N j F j N+l = 1. (135) 
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Thus, by defining = to mean equality up though and including terms of order t we see 
immediately from the form ([T2]i with (|42j! for Kq } and (|46]l, (@7) and (g9]l for the with 
<n < N that the terms though order in f^ N {t) are 

which vanishes by use of J135I ). This derivation has made no use of c^ N . This term contributes 
only to order t N+l and may be determined from the normalization amplitude (12TV This 
provides an alternative to the derivation of (l48~t of Appendix B. 

To prove the cancellation for J^l,(t) we note that because of the term C„ (N; t) in 

C%\N; t) for m = 0, 1 , 2 in (|53) we may use the expression (UJ) and d42l-(l49l for /^(f) 
in the form 

f^C ( i\N;t)-F 2 N m F™ +l = (137) 

m=0 

Thus from ([T3K (fTOt and (1137b we obtain an alternative form for fS'Lff) of 

O0= (| + ^0'>VO'> +^ 2 -i^;0-4-^ +I , (138) 

^ J m=0 

where 

2 /«\ 2N+\-m 

Ci(N; t) = (-1)"+ 1 • - ■ " • ^ ■ J] c®, (139) 

We have already demonstrated by use of (11361 > that ffjfL(t) vanishes though order f^. Therefore 
using the expressions (l55ll-(l5"8l for cj 3 .', which are all valid through (at least) order t N and the 
definition® of /^(f) we find 

which vanishes by use of the Wronskian relation (1135b with n — 3. 

We have thus demonstrated that f^ N (t) I t N ^ 2 vanishes to order t N as t — > 0. However we 

see, from the orig inal integral <TT8>, that in fact jfUt)/^ 12 is of order t N+2 . Therefore the 
coefficient of t N+1 must also vanish. This is not proven by (11401 i. However the coefficient Cj 3 ^ 
has not been used in the derivation of ( 1140b and the choice of Cj 3 ^ to make the coefficient of 
t N+l vanish provides an alternative derivation of (l59"V 

6. Factorization for f£> with n > 4 

In principle the methods of differential algebra of the previous sections can be extended to 
form factors fuL(f) with n > 4. However, the complexity of the calculations rapidly increases. 

For fxx(t) there are In + 1 polynomials C® n \N; t) and since from [ 15 1 we find that for 
N > \ 

L 2n+l ■■■ Li-L^-f^Jit) = 0, (141) 



~u 2 (N+ 1) 



2-7 



■u 2 (NyF 2 - J Fj._ 



(136) 



« 2 (JV + 1) 



3-7 



«2(JV)' • F 



3-7 



F ; 1 

1 N+l I 



(140) 
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where Lf. is a linear differential operator of order k, the polynomials C^"\N; f) will satisfy a 
system of 2n + 1 coupled differential equations where the maximum derivative order is n 2 . 
These equations can be decoupled into 2n + 1 Fuchsian ODEs which genetically have order 
n 2 (2« + l). 

Similarly for f^V~ l \t) we found in [ 15 1 that 

L 2n+2 ■■■ L 4 -L 2 -f£" N +r, (t) = 0, (142) 

and thus the 2n + 2 polynomials C%" +r> (N; t) satisfy inhomogeneous coupled equations of 
maximum differential order n(n + 1) which for N > 1 are genetically decoupled into Fuchsian 
equations of order 2n(n + l) 2 . 

We have obtained for f^ N (t) the 20-th order ODEs for C ( *\N; t) in the cases N - 
1, ■ • • ,10 and will illustrate the new features which arise by considering the case m — 4. 

We find by use of Maple that (at least for low values of N) the operator Q.f\N; t) has a 
direct sum decomposition 

Glf{N\t) = Mf\N)®M^(N)®Mf^(N)@Mf\N), (143) 

where M ( ^ n (N) is order k and is homomorphic to the symmetric k - 1 power of #2 (AO 

Mf{N)-jf{N\t) = G%\N; t) ■ Sym 6 {0 2 {N)), (144) 

M^(N) ■ jf\N;t) = Gf\N; t) ■ Sym 4 (<9 2 (A0), (145) 

m£J(A0 = Sym 4 (O 2 (A0), (146) 

Mf{N)-J a {N;t) = G ( *\N;t)-Sym 2 (0 2 (N)), (147) 

where the intertwinners /^ 4) (2; f) and Gj, 4) (2; t) are linear differential operator of order m. The 
intertwinners J„ (2; t) in (11441 i- (11471 i. are explicitly given in Appendix F Further examples of 



intertwinners are given in |Appendix F| These differential algebra exact results (in particular 
(!1441 l-( fT47T l) are the illustration of the canonical link between the palindromic polynomials 
and the theory of elliptic curves. 

Direct sum decomposition^ have been obtained for £l„\N; t), Q$(N; t) and O^CA 7 ; t) 
and we conjecture that this occurs generically for all Q„ \N; t). Taking into account the 
homomorphism of 2 (N; t) and 2 (N + T, f), and recalling, for instance, subsections 14. 2t and 
14.31 it may be easier to write direct sum decomposition formulae in terms of sum of symmetric 
products of 2 (N; t) and 2 (N+ 1 ; t). In order to extend these results, beyond these few special 
cases of q!£\N; t), a deeper and systematic study of the homomorphisms is still required. 

From an analytical viewpoint, a complication which needs to be understood is how to use 
the solutions of the homogeneous operators Q. ( ^\N; t) to obtain the polynomial solution of the 
inhomogeneous equations. The first difficulty here is that for C„(N; t) the inhomogeneous 
terms are large polynomials, of order 100 and higher. Moreover, the orders of palindromy 
point of the cf{N; t) with N = 1, • • • , 10 are all larger than the order f A ' +4 where the solutions 
of the homogeneous operators M^(N; t) have their first logarithmic singularity. Consequently 
linear combinations of solutions must be made which cancel these logarithmic singularities 
at t N+4 to give sets of solutions to Q.f\N; t) which are analytic up to the order of the first 
inhomogeneous terms. Thus the determination of the correct linear combination of solutions 

t Note that in direct sum decomposition like 11431 . some ambiguity may occur with terms like M^(N) ffi A/j j(iV) 
where M» ?(jV) and Mil(N) are both homomorphic to a same operator (here Sym 4 (02(^V)))- 



23 



of the operators M^(N;t) is significantly more complex than was the case for C„\N;f). 

Exact results for the C* 4) 's, based on the Wronskian cancellation method of Sec. [5] and valid 
for any value of N are displayed in Appendix G These are exact results for the palindromic 
polynomials in terms of Fn and U2(N), namely two hypergeometric functions associated 
with elliptic curves. Thus, these analytical results can also be seen as an illustration of the 
canonical link between our palindromic polynomials and the theory of elliptic curves. They 
confirm the deep relation we find, algebraically and analytically, on these structures with the 
theory of elliptic curves. In a forthcoming publication we will show that the relation is in fact, 
more specifically, a close relation with modular forms. 



7. Conclusions 

In this paper we have proven the factorization, for all N, of the diagonal form factor f]y N {t) 
for n = 2, 3 previously seen in [ 15 1 for N < 4 and provided a conjecture for n = 4. Besides 
new results like the quadratic recursion (|5H , or non trivial quartic identities (like (|55ll- d58Tl), 
one of the main result of the paper is the fact that, introducing the selected hypergeometric 
functions F N , which are also elliptic functions, and are simply related to the (simplest) form 
factor k, jL the form factors actually become polynomials of these F^s with palindromic 
polynomial coefficients. The complexity of the form factors, is, thus, reduced to some 
encoding in terms of palindromic polynomials. As a consequence, understanding the form 
factors amounts to describing and understanding an infinite set of palindromic polynomials, 
canonically associated with elliptic curves. 

We also observe that all of these palindromic polynomials are built from the solutions of 
the operator C>2(N), and, therefore, are all properties of the basic elliptic curve which underlies 
all computations of the Ising model. There is a deep structure here which needs to be greatly 
developed. The differential algebra approach of the linear differential operators associated 
with these palindromic polynomials is found to be a surprisingly rich structure canonically 
associated with elliptic curves. In a forthcoming publication, we will show that such rich 
structures are closely related to modular forms. 

Analytically, the conjecture and the Wronskian method of logarithm cancellation can 
be extended to large values of n, but the method of proof by differential equations becomes 
prohibitively cumbersome for n > 4. This is very similar to the situation which occurred 
for the factorization of correlations in the XXZ model where the factorizations of 1 19|-[25| 
done for small values of the separation of the spins by means of explicit computations on 
integrals was proven for all separations in ll30l by means of the qKZ equation satisfied by 
the correlations and not by the explicit integrals which are the solution of this equation. This 
suggests that our palindromic polynomials may profitably be considered as a specialization 
of polynomials of n variables. Moreover, if the two conjectures presented in the introduction 
are indeed correct, then such kind of structures could also have relevance to the 8 vertex 
model and to the higher genus curves which arise in the chiral Potts model. Consequently the 
computations presented here could be a special case of a much larger modularity phenomenon. 
This could presumably generalize the relations which the Ising model has with modular forms 
and Calabi-Yau structures ll3"T| . 
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Appendix A. Form factors in the basis F N and F, 



By use of the contiguous relations for hypergeometric functions the examples given in lfl5l of 
ff? N (f) expressed in terms of the elliptic integrals K(t l l 2 ) and E(t 1 ^ 2 ) may be re-expressed in 
terms of the functions Fn and F^+i- Several examples are as follows 



/ao = j-Fo-Fu (A.1) 
J$ = \ -\(t + l){2t 2 + t + 2).F\ 



3 2 3 4 
+ s-t- (At 2 + 5 t + 4) • F x ■ F 2 - -7 t 2 (t + 1) • F 2 (A.2) 

2 5 2 7 



jfl = 1 -^(r+l)(64r 4 +16f 3 +99f 2 + 16f + 64)-F 2 



5 2 

t ■ (64f 4 + 88 r 3 + 105 1 2 + 88 1 + 64) ■ F 2 ■ F 3 



2 s -3 
5 4 
2 7 -3 2 



•f 2 -(f+ l)(2f 2 + f + 2)-F^, (A.3) 



H2) 



3 



f 

J 3.3 2 



• (f + 1) (576 r 6 + 96 r 5 + 730 f 4 + 425 f 3 + 730 f 2 + 96 1 + 576) • F 2 

7 2 



2 I2 -3 

74 



• t (768f 6 + 928 f 5 + 1240 f 4 + 1455 f 3 + 1240 f 2 + 928 1 + 768) • F 3 F 4 



7 4 

^— ^ • t 2 (t + 1) (64 1 4 + 16 f 3 + 99 1 2 + 16 f + 64) • F 2 . (A.4) 



For with TV = 0, ■ • • ,4 



/o,o = 23 ~— !l ^->V I--/-/V-F,- <A.5» 

< 3 > „ /-CD 



1 



/■I j i 
■M.l 



3 



f 1 / 2 3 f'/ 2 2 3 -3 

+ (8f 2 + 15 t + 8) • Fx 1 ■ F 2 - ^rt 2 (t + l)FxF 2 + If 3 . ^ (A 6) 

f (3) „ f (D 

J2,2 _ 1 -'2,2 

~T ~ 2~3 ' ~T 

- (1 + f) (2 6 • 3 • 7 f 4 + 1 136 f 3 + 3229 ? 2 + 1 136 1 + 1344) • F 2 3 

5 2 

(2 5 ■ 3 2 ? 4 + 596 f 3 + 859 1 2 + 596 t + 2 5 • 3 2 ) ■ F 2 2 • F 3 



2 11 -3 



o ;o 32 (r + l)(3r 2 + At + 3) f 2 • F 2 ■ F 2 + ^T^ 1 ' (3f 2 + 8 r + 3) • F 3 \ (A.7) 
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f (3> , f m 

^3,3 _ J J_X3_ 



- 2U l (t + 1)(2 7 • 3 3 • 5 2 f 6 + 49680 f 5 + 153306 t 4 + 160427 f 3 

+ 153306 f 2 + 49680 1 + 2 7 • 3 3 • 5 2 ) • f 3 
7 2 

+ -rr — rf(2 10 ■ 3 2 ■ 5f 6 + 79200 r 5 + 128104 t 4 + 168593 f 3 
2 ie . 33 

+ 128104 f 2 + 79200 f + 2 10 • 3 2 • 5) • F\ ■ F 4 

r-,A 

r(t + l)f 2 (2 4 • 3 2 • 5t 4 + 670 f 3 + 1763 f 2 + 670 f + 2 4 • 3 2 • 5) • F 3 • F\ 



■2 



7 6 

221 3 4 ? 3 (2 6 • 5f 4 + 740 f 3 + 1407 f 2 + 740 f + 2 6 • 5) • F\, (A.8) 



f (3) f (D 
J44 ij J44 



- 2^-^ (f + 1) (2 14 • 5 • 7 • 13f 8 + 3254272 1 1 + 1 1474624 f 6 + 8672032 f 5 

+ 20423231 f 4 + 8672032 f 3 + 1 1474624 f 2 + 3254272 f + 2 14 • 5 • 7 • 13) • f\ 

33 

+ — f(2 12 • 3 • 5 2 • 7r 8 + 3334912 1 1 + 4845120 1 6 + 7068720 f 5 

2 23 • 5 

+ 8865649 t 4 + 7068720 f 3 + 4845120 f 2 + 3334912 1 + 2 12 • 3 • 5 2 • 7) • F\ ■ F 5 

3 6 



2 20 . 5 3 

)37 

3 8 



f 2 (f + 1)(2 4 • 3 2 • 5 • 7 2 f 6 + 26292 f 5 + 69377 1 4 + 78580 f 3 
2 + 26292 1 + 2 4 • 3 2 • 5 • 7 2 ) • F 4 ■ Fj 
t 3 (2 3 ■ 3 3 • 5 • 7 f 6 + 2 3 • 3 3 • 7 • 1 If 5 + 28413 f 4 + 46432 f 3 



2 20 . 5 3 

+ 28413 f 2 + 2 3 • 3 3 • 7 • 1 1 1 + 2 3 • 3 3 • 5 • 7) • F\. (A.9) 

The coefficients which are not given in factored form all contain large prime factors. 
For with N =0,1,2,3 

C = rffi-w^-'-n ^-fFi-Fl (A.10) 

- ^t(t + l)(8f 4 + 18f 3 + 35f 2 + 18r + 8) • Fl ■ F 2 

34 

+ — • f 2 (8f 4 + 28f 3 + 45f 2 + 28f + 8) • F\ ■ F\ 

-|l-f 3 -(f+l)-(4f 2 +llf + 4)- J F 1 • J F 3 + |lf 4 (f 2 +4f+l)- J F 4 , (All) 
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f< 4 > - _I 4. — f( 2 ) 
2,2 ~~ 3 3' 2 ' 2 



+ (2 14 f 10 + 40960 t 9 + 84480 f 8 + 136640 1 1 + 176180 1 6 

+201075 f 5 + 176180 f 4 + 136640 f 3 + 84480 t 2 + 40960 1 + 2 14 ) • F 4 

5 2 

- 214 32 f (t + 1) (2 13 f 8 + 13312 f + 29504 f 6 + 36320 f 5 + 45337 t 4 

+36320 f 3 + 29504 t 2 + 13312 f + 2 13 ) • F\ ■ F 3 
5 4 

+ — - — - t 2 (2 12 f 8 + 1 1264 f + 21760 f 6 + 31576 f 5 + 36209 f 4 
2 17 • 3 2 v 

+31576 f 3 + 21760 f 2 + 1 1264 1 + 2 12 ) ^1^1 
5 6 



■3 4 



f 3 (f + 1) (2 8 f 6 + 480 f 5 + 906 t 4 + 979 f 3 + 906 f 2 



+480 1 + 2 8 ) • F 2 • 

5 8 
2 i5 . 34 



58 

+ f 4 (2 5 f 6 + 96 f 5 + 177 1 4 + 224 f 3 + 177 f 2 + 96 1 + 2 s ) • F 4 , (A. 12) 



f (4) _ _ 'i . JJ_ f (2) 

^3,3 - 2 3 2 • 3 ' 3,3 



+ ,} - (2 13 • 3 4 • 7f 14 + 10838016 f 13 + 19643904 f 12 + 34169856 t u 
2i6 . 34 V 

+50403584 f 10 + 62791680 1 9 + 73309425 f 8 + 79935700 1 1 + 73309425 f 6 



+62791680 f + 50403584 f + 34169856 f + 19643904 t L + 10838016 1 
+2 13 -3 4 -7)- J F 4 
7 2 

- 2l9 34 t(t+\) (2 14 • 3 3 • 7f 12 + 4257792 f 1 1 + 9547776 t w + 13813120 1 9 
+ 1 934 1 1 20 f 8 + 2 1 399090 1 1 + 24976435 1 6 + 2 1 399090 1 5 + 1 934 1 1 20 1 4 

+ 13813120 f 3 + 9547776 1 2 + 4257792 1 + 2 14 • 3 3 • 7) • F\ ■ F 4 
7 4 

+ -r, — ? f 2 • (2 15 • 3 2 • 7 t 12 + 4988928 f 11 + 9680384 f 10 + 15992320 1 9 

2 25 . 3 5 V 

+2 1 863 1 20 f 8 + 26325960 1 1 + 285270 1 5 1 6 + 26325960 1 5 + 2 1 863 1 20 1 4 
+ 15992320 f 3 + 9680384 1 2 + 4988928 f + 2 15 • 3 2 • 7) • F\ ■ F\ 

7 6 

- 22? 34 f 3 • (f + 1) (2 14 • 3 • 7f 10 + 501760 1 9 + 1 191680 f 8 + 1548640 f 7 

+2065400 f 6 + 2169745 f 5 + 2065400 1 4 + 1548640 f 3 
+ 1 191680 1 2 + 501760 f + 2 14 • 3 • 7) • F 3 • F 3 4 

7 8 

+ --: — j t 4 ■ (2 12 • lt m + 71680 1 9 + 147840 f 8 + 235040 1 1 + 299555 1 6 
2 3i . 3 4 V 

+339180 f 5 + 299555 1 4 + 235040 f 3 + 147840 1 2 + 71680 f + 2 12 • 7) • F 4 . (A. 13) 
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For with N = 1,2,3 

AS) o2 f (D A3) 

J\,\ _ J\,\ J\,\ 

7^ ~ J'7^ + 7^ 

+ (t + 1) 2 (2 6 + 136f 3 + 159f 2 + 136? + 2 6 ) • F? 

2 6 • 3 • 5 1 

- ^r(r + l)(2 5 f 4 + 80f 3 + 99t 2 + 80r + 2 s ) • F 4 • F 2 

3 3 

+ ^2f 2 (2 7 f 4 + 368f 3 + 483f 2 + 368f + 2 7 ) • F 3 • F\ 

- ^t\t +l)(4f 2 +5f + 4) ^1^1 

+ |^f 4 (8f 2 + 13f + 8) • Fi ■ F\ - J^(f + 1) • f 5 • F\, (A.14) 



fl,! _ 137 f 2 J 3 ^2,2 
~T ~ ~2 3 -5 ' T~ + 2 ' ~T 

+ 218 1 - (8 f 2 + 7 1 + 8) (2 9 • 3 • 61 f 8 + 241856 1 1 + 508200 f 6 + 708609 f 5 
+780244 f 4 + 708609 f 3 + 508200 f 2 + 241856 f + 2 9 • 3 • 6l) • F 5 2 

- 218 32 f (f + 1) (92160 f 8 + 239360 f 7 + 540576 f 6 + 723924 f 5 
+868861 f 4 + 723924 f 3 + 540576 1 1 + 239360 1 + 92160) • F 4 • F 3 



5 4 

1 (90624 f 8 + 338816 f 7 + 743304 t h + 1 122432 f 5 + 1278697 f 4 



+ f 2 

2 20 • 3 3 

+ 1 122432 f 3 + 743304 f 2 + 338816 f + 90624) • F 2 3 F 3 2 
5 6 

- 2l? 34 f 3 (f + l)(l392 f 6 + 4010 f 5 + 6983 f 4 + 8136 f 3 

+6983 f 2 + 4010 1 + 1392) • F 2 • F 3 
5 8 

+ — — - f 4 (684 f 6 + 2752 f 5 +5161 f 4 + 6240 f 3 + 5161 1 2 + 2752 r 
2 20 . 3 5 V 

+684) • F 2 • F 4 

5 9 



2 19 . 36 



f 5 + 1) (42 f 4 + 133 f 3 + 167 1 1 + 133 f + 42) • F\, (A. 15) 
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AS) AD A3) 

J 3,3 _ i^-' ^3,3 ^ ■'3,3 

+ — — ; ; — (2 16 • 3 4 • 5 • 17 f 14 + 1377976320 1 13 + 3016452096 f 12 
2 20 . 3 s . 5 I 

+5930641920 f 11 +9308313280 f 10 + 12328157240 f 9 + 14834544515 f 8 
+ 15849843292 r 7 + 14834544515 f 6 + 12328157240 f 5 + 9308313280 1 4 



+5930641920 f 3 + 3016452096 t 2 + 1377976320 1 + 2 16 ■ 3 4 ■ 5 • 17) • F\ 

- 222 35 f(f + 1)(2 19 -3 3 -5t n + 15151 1040 f 11 +351000576 t w 

+605214208 1 9 + 835692208 f 8 + 1025976166 f + 1 1 12168875 1 6 

+ 1025976166 f 5 + 835692208 f 4 + 605214208 f 3 + 351000576 t 2 

+ 151511040f + 2 19 -3 3 -5)-F 4 -F 4 
74 

+ -= — - f 2 (2 18 • 3 3 • 5 2 f 12 + 572129280 f 11 + 1334317056 t m 
2 29 . 35 V 

+2446757888 f 9 + 3545541888 r 8 + 4425343776 f 7 + 4784608975 1 6 
+4425343776 1 5 + 3545541888 r 4 + 2446757888 r 3 + 1334317056^ 

+572129280 1 + 2 18 • 3 3 ■ 5 2 ) ■ F\ ■ F\ 

76 

- — — - r 3 (r + 1) (2 13 • 3 • 5 • 7 1 10 + 2007040 f 9 + 4885888 r 8 + 7228048 f 

227 . ^5 v 7 \ 

+9666130 ? 6 + 10423545 f 5 + 9666130 f 4 + 7228048 ? 3 + 4885888 f 2 

+2007040 f + 2 13 • 3 • 5 • 7) • F] ■ F 3 4 

78 

+ -r A — r t 4 (2 14 • 5 • 13 t w + 3665920 f 9 + 9078784 f 8 + 15185664 1 1 

2 34 . 35 V 

+20375540 1 6 + 22605185 f 5 + 20375540 1 4 + 15185664 f 3 + 9078784 f 2 

+3665920 r+2 14 -5- 13) -/V^ 4 . 
710 

-— — ; — f 5 (f + 1)(2 13 -5f 8 + 104960 1 1 + 2671 36 1 6 + 3 19904 1 5 

2 35 . 3 5 . 5 '\ 

+436441 1 4 + 319904 f 3 + 267136 1 2 + 104960 f + 2 13 • 5) • F|. (A. 16) 



Appendix B. Polynomial solution calculations for cf\N; t) 

(2) 

We here give explicitly the calculational details for C\ (N; t). 

Using the form (l43l in the inhomogeneous equation (l83l we find the recursion relation 
for the coefficients c (2) n for n + N, N + 1, N + 2, N + 3 

2n-(n-N)(n-2N - I) ■ cf n 

- {2n 3 - 6JVn 2 - 2(4 + JV - 2JV 2 )n + 5 + 6N} ■ c®_j 

- {2n 3 - 6(3 + JV)n 2 + (46 + 3AN + AN 2 )n - 35 - 38JV - 8JV 2 } • c®_ 2 

+ 2 (n - 2)(n -2N - 3)(n - iV - 3) • c®_ 3 = 0. (B.l) 
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where by definition c\. n — for n < — 1. This recursion relation has four terms instead of the 
three terms in the corresponding relation (|92l for c? . We note that, if we send n — > 2N-n+3 
in ( IB. lb . we see that c^ 2 ? and c®. satisfy the same equation. Since the coefficient of Cj 
vanishes for n = 0, the term cfi is not determined from (IB.lt and by convention we set 

c (2) =1 

L l;0 1 

Following the procedure used for cf\N; t) we note that equation d84l i will be satisfied 
to order t* 1 if we choose the Cj for < n < A/ - 1 to be the corresponding coefficients in 
r 1 ■ u 2 (N; t) ■ u 2 (N + 1; t) and hence (|47j follows. 

The inhomogeneous recursion relations for n — N, N + 1 are 

A (2) {-(2N 2 + 2N-5)- cf^_j + (8A/ 2 + 8A/ - 35) • c (2 ^_ 2 

-6(N-2)(N + 3)-cf^_ 3 ] = -2N 2 (2N+\) 2 A 2 N , (B.2) 



A (2) {-2 N(N+l)- c® +1 -(2N+1) 2 - c% 



+ (6N 2 + 6N-5)- c® _, + 4 (tf + 2)(tf - 1) ■ c®_ 2 } 

(2N + l) 2 (4A/ 3 + 4N 2 - 4N - 1) A ; 



w (B.3) 



4(AT+1) 

and the relations for N + 2, N + 3 are identical with iV, AT + 1, respectively, with the 
(palindromic) replacement 

c <2) — » c (2) (B4) 
If there were no inhomogeneous term ( IB. 2) would be a new constraint in the coefficients 

(2) 

c. for n — N - 1 , N — 2, N - 3. However this constraint does not hold (because the solution 
to the homogeneous equation has a term f A ' +1 In f)- 

The normalizing constant A (2) can be evaluated from (IB.2l i and the sum on the LHS of 
( IB. 2b is evaluated the same way the corresponding sum was for C 2 2 \N; t), by comparing with 
the full solution r l U2(N; f) u 2 (N + 1; t) of the homogeneous equation. Thus we find 

- (2A/ 2 + 2AT -5) • cf^CA 7 ) + (8A/ 2 + SN - 35) • c® _ 2 (A0 

-6 (A/ + 3) (TV -2)- c®_ 3 (A0 = -2 JV 2 (JV + 1) 4, (B.5) 

and, hence, we find from dB.2| > 

A (2) = JV0y. (B.6) 

It remains to compute c {2) N from dB.3l l. We obtain the palindromic solution by requiring 
that c (2 ^ +1 = cyL_j and thus (IB .3b reduces to 

\ 2 . -( 2 ) ^/i at 2 /i a? _ s\ . J- 2 ) /lew _i_ ww _ n . J® 



(2N + iy ■ + (4JV Z + 4JV - 5) • e^_, + 4(JV + 2)(iV - 1) • c^_ 2 

. A AT _ 1 ^ 



(2W + l) 2 (4A/ 3 + 4A/ 2 - 4JV - 1) /I 2 



4(A/+ 1) 



(B.7) 



An equivalent and more efficient method for evaluating c (2 ^, which avoids the need to evaluate 
the sums on the LHS of ( IB.7I ), is to directly evaluate C (2) (N; t) in terms of Cf\N; t) by use of 
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the coupled equation (|82] i. From this we find 

= 4 +Y j ak(N)-a N - l - k (N), (B.8) 



N-l 

jn ( 

k=0 



and, by explicitly evaluating the sum in ( IB. 81 ). we obtain the result d48l . Finally, the Cj 2) n for 
iV + 1 < « < 2./V are determined from the palindromy of (IB.lt . 



Appendix C. Coupled differential equations for Cn\N; t) 
The four coupled differential equations for C$(N; t) are 



2Af + 1 (3) (N+l)(2tN+l+t) (3) 

2-(f-l)f W ? 2 -(2JV +!)(?-!) 1 1 



(jV+l) 2 (2fjV + 3) (3) Q (jV+l) 3 (fjV-f + 3) (3) 
f 2 (2A^+ l) 2 (f- 1) 2 W (f-l)(2JV+l) 3 f 2 3 ( 



tN + 2t-N-\ ± c o) ( . 2 (N+l)(tN-N + t) d_ (3) 
(f- l)r A U f(f- 1)(2JV + 1) * 1 



(A^m^ . l c (3) (f) + 8 (N + l?( t N-N + 2- t) ± 
t(t- l)(2N+ l) 2 * 2 f(2JV+ l) 3 (f- 1) dt 3 

* A „ *+l _ ^ . (jV + l) 2 _ J 2 



j cr '(t) + 2 ■ Cr '(t) + 4 — — ■ C { >(t) 

+ dfi {> 2N+1 dt 2 1 (2N + l) 2 dt 2 2 

+ 8 (jV+1)3 4ld 3) (f) = ^ r^ 1 -(2iV + l). B (N), (C.l) 



(jV+l)(2jV + 6fjV + 3f + 2) (3) 

f 2 (2jV+l) 2 1 W 

o (tf + l) 2 (4 riV + 4 + 5 + _ c(3)(f) 



t 2 (2N+ l) 3 

l (jV+1)3(6jV + 2 f 2 +9 - 2f) .cf(Q + 6^cf( 
(2JV + l) 4 f 2 3 (if 

(jV+l)(5fjV + jV + 3f+l) ^ c (3) f . 

t(2N + l) 2 A 1 



o (jV+l) 2 (2jV + 4fjV + f + 4) d » 
+ 8 f (2JV+l) 3 ^ 2 ( ° 

f(2JV+ l) 4 <ft 3 (2JV+ l) 2 dt 2 1 W 

(;v + i) 2 (;-i) rf 2 (3) (jv + i) 3 (f-i) d 2 (3) 

+ 16 (2JV+1) 3 (0+48 (2JV+1) 4 'a 2 ^ 3 ( ° 

= - t N ~ l ■(2N 2 t + t + 4-tN-2N-2N 2 )-B Q (N), (C.2) 
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1 V 7 f(2JV+ l) 2 2 (2JV+ l) 2 c/f 2 W 



+ 4 



(AT + l) 2 (12 AT 2 -l6tN - 2t N + 30 N + 2 r + 18 - 17 f) 



(2N + Iff 
(t-l)(N+l) 2 (tN + 5N-t + 5) d 



3 (2AT + l) 2 

4 (JV+ 1) 



t(2N+ l) 4 



(f- l) 2 (jV+ l) 2 
(2JV+ l) 4 



— C (3) (f) 



• f^ +1 • B (N), 



(C3) 



(3) (iV + l)(4iV + 2riV + 4-0 (3) 

W t(2N+l) 2 1 W 

(N +1) 2 (4N 2 + 8N 2 t-10t 2 N +10tN +12N-t-4t 2 + 8) 



Cf{t) 

(N+ 1) 3 (4JV Z -2t z N- 10 tN+ \6N + 2t l + 13 - 14 1) ^ (3) 



(2N + l) 4 f 2 



+ 48 — — • C ^(?) 



(2JV + l) 5 f 2 
(r-l)(AT + l) £ (3) 
(2JV+1) 2 & 1 W 



(t - \)(N + X) 2 (5tN + 3N + 2t + 3) d 

+ 16 . ,, 4 tA a) 



<3), 

t(2N+l) 4 dt 



(jV+l)3( f -l)( f jy + 3jV + 4-0 4 (3) 

+ 96 ■ — '(t) 

t{2N+\f dt 3 



(2JV+1) 4 dt 2 2 ' (2JV+1) 5 t/f 2 3 

= 3 t N -(3N(t-l) +2t-l)-B (N), (C4) 



where Bq(N) is given by (II 1 lb . 

Appendix D. The ODE and recursion relation for cf\N; t) 

The ODE for cf(N; f) can be found by carefully using the four coupled ODEs (lCTt-(|C4li. 
First use (IC.3b to solve for Cf'(N; t) and then use this is in Equation ( IC.4b in order to solve 
for Cf(N; t). Next, use both Cf{N; t) and Cf\n; t) in ED and to produce ODEs of 
orders four in cf (N, t) and five in cf\N; t) in dC.U and orders three in C? (N; t) and four in 
Cf (N; t) in Oil- 
In the new ( IC.U . the fourth derivative of C^iN; t) can be solved in terms of the other 
derivatives, and likewise in the new (lC.21 i. the third derivative of Cf (N; t) can be solved in 
terms of the other derivatives. Taking the derivative of the expression for the third derivative 

(3) (3) 

of C 2 (N; t) and equating it to the expression for the fourth derivative of C\ (N; t) we find an 
alternate expression for the third derivative of C\ (N; t). Finally, equating the two expressions 
for the third derivative of (N; t), a full cancellation of all of the derivatives of C^{N; t) 
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takes place, leaving a fifth order ODE in terms of only Cf\t) 

4 [2 (N - l)(2N + l)(3N + l)(N + l)f 4 

-(2N + 3)(36N 3 -7N 2 -69N-32) f 3 
+4 (N + 2) (36 N 3 - 10 N 2 - 1 16 N - 69) t 2 
-(2 N + 5) (60 N 3 -23N 2 - 275 AT - 188)f 
+ 18(2^ + 3)(A? + 3)(^-3)(A^+ 1)] •Cf(f) 
-8[(AT-1)(2AT+1)(3AT+ 1)(AT + l)f 4 

-(-130N + 40Ar 3 -47 + 24 N 4 - 73 N 2 ) r 3 
+2 (l8N 4 - 129 +45 A? 3 - 1 1 3 N 2 - 270 n) t 2 
-(-253 N 2 - 422 + 24 N 4 - 740 N + 80 N 3 ) t 

+ (N+l)(6N 3 + 19N 2 - 114N-211)]-* - — cf (t) 

+ 20(f- l)[2 N(N- l)(N+ l)-f 3 

-3 (2 A? 3 - 4 A? 2 - 8 N - 3) ■ t 2 + 3 (-13 + 2 JV 3 - 8 N 2 - 24 AO f] 

d 2 



-2 (N -9)(N + 2) (N + 1)] • t ■ —C^'{t) 



,2 u r (3), 
dt 2 

j3 

+ 40(f- 1) 2 [(N- l) 2 r 2 - (4-N+ 1 +2N 2 )t + (N + 5)(N + 1)] ■ t 3 ■ ^ c rW 



j4 j5 

- 40 (f - l) 3 [(AT - l)f - A? - 1] f 4 ■ — rd 3 V) + 8 (f - l) 4 • f 5 • — C' 3) (f) 

ar dp 

From this differential equation we obtain the recursion relation for the coefficients c\. n 

and the normalization constant A ( p defined by the form (|53T >. where by definition c?^ = for 
n < -1 

Af ■ {8n(2N-n)(N-n)(N + n)(3N-n) ■ cf n 
+ 4 (2N +1 -2n){2-l n + 1 N - N 2 + 4« 4 - 12 N 3 - 8 n 3 + 24N 3 n 
-4 A? 2 n + 1 1 « 2 + 4 AfV - 16 Nn 3 + 24 Nn 2 - 22 Nn) • c i3) , 

- 16 (N + 1 - n) (9 - 22 n + 22 N + N 2 + 3 n A - 1 8 N 3 - 12 n 3 + 1 8 N 3 n 
-6 N 2 n + 23 n 2 + 3 N 2 n 2 - 12 Nn 3 + 36 Nn 2 - 46 Nn) • c {3) . 

) 3;«-2 

+ 4 (2 A? + 3 -2n)(32-69n + 69N + 7 A? 2 + 4n 4 -36N 3 -24n 3 + 2AN 3 n 
-\2N 2 n +59« 2 +4Af 2 « 2 -16A^n 3 + 72Afa 2 - 118 Nn) • c (3) , 

/ 5\n— j 

- 8 (n - 2)(2N + 2 - n)(N + 2 - n)(N - 2 + n)(3 N + 2 - n) ■ c (3) 

3(2N + l) 6 

- (8 n ,N ~ ^n,N+2) ■ —777 777" ' Bq. (D.2) 

(A 7 + l) 3 

We note by sending n — > 2A 7 - n + 2 that and Cj 3 l v _ 2 _ satisfy the same equation. 



3 

(3; 
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For n — ( ID. 21 ) is identically zero for any c?g which we set equal to unity by convention 

,(3) 
"3;n 

>geneous equation 

1WC111 Ul C 3;A 

term, the coefficients c?^ for n = N - 4, N-3,N-2N-1 would have to satisfy a 
non trivial constraint. This constraint does not, in fact, hold and is the reason that the 
homogeneous equation has a term f A ' +3 In t. However, with a nonvanishing inhomogeneous 
term, the equation for n — N determines the normalization constant. 
For n = N + 1 the equation (ID. 21 ) reduces to 

8 (N + 1) (N - 1) (2 N + 1) (2N - 1) • (c® +1 - cfj,_ 4 ) 
-8 (2^V+l)(2^V-l)(2^V 2 -l)-(4^-4^_ 2 ) = (D.3) 

which will be satisfied by the palindromic property 

c ( ? = 4 3 > , (D.4) 

with n — N + 1 and n - N. Finally, the c\. for N < n < 2N - 2 are determined from the 
palindromy of dD.2t . 

Appendix E. Homomorphisms for C„ 3) (W; f) and cf\N; t) 

The fifth order operator Mt'(N; t) in the direct sum decomposition (1122b of Q^{N; t) has the 
homomorphism (in terms of the operator LziN)) 

M^\N)-tf\N;t) = G^(N;t)-Sym 4 (L 2 (N+l)), (E.l) 
where the intertwinners / () (N; t) and G (N; t) are: 

jO\ N ; t) = f ^. (f _l). f .( A -l^) 

= f +l -((f- l)-r-D, -(2 AT + 2(N + 1)) ), (E.2) 
Gf^O = ,™.(,-i). f .(z>,-£^), (E.3) 



For < n < N- 1 the rhs of (1D.2I) vanishes and hence the cl 3) are identical with the coefficients 



( 11 18t of the solution (11 16b ot the homogeneous equation. 

For n - N the coefficient of CjjL vanishes, and thus if there were no inhomogeneous 



where 



4 = (f _ 1) 2(2 W + l). r 2( i V + l) ; (E4) 

nfi _ (r + l)( f -l) 4;v - 3 

- jWTb Pn, (E.5) 

= (4N + 3)-(3N + 2)- (f 2 + 1) + 2(20N 2 + 15 N + 2) • f 

= (4JV + 3)-(3JV + 2)-(f+ l) 2 (E.6) 



(3) 

The homomorphism for Mi" (N; t) is 



+ 4 (2(2 N + 1) (AT - 1) +A0 • t. 

2 



M?\N; t) ■ £\N; t) = Gf (N; t) ■ Sym 4 (L 2 (N + 1)), (E.7) 
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where the intertwinners jf\N; t) and G^(N; t) are 



jf\N;f) = f +2 •((*- l)-f D, + 2 (N+ l)-t + 2N), (E.8) 

d\n(R B J-(N+ l))) x 
Gf\N; f) = f" +2 • (f - 1) • f ■ A - CE.9) 



where R?, is exactly the in (IE.5b - 



Appendix F. Homomorphisms for \N; t) 

Many exact results have been obtained on the intertwinners occurring in (I144K (1145K (1146K 
dl47| > . Let us display the simplest ones. 
For Jq\N; t) we have 

/< 4) (2;f) = f 2 -(f+l)-(2f 2 + f + 2), 

J (4) (3;t) = f 2 -(f+l)-(64f 4 + 16f 3 +99r 2 + 16f + 64), 

y< 4) (4; f) = t 2 ■ (t + 1) ■ (576 f 6 + 96 f 5 + 730 t 4 + 425 f 3 + 730 f 2 + 96 1 + 516) , 

f 4) {5\t) = t 2 ■ (t + 1)- (16384 f 8 + 2048 1 1 + 19264 t 6 

+6608 f 5 + 28861 f 4 + 6608 f 3 + 19264 f 2 + 2048 r+ 16384). (F. 1) 

For j[\N; t) we have 

2 f 3 • yj 4) (2; t) = (f - 1) • 7 ( (4) (2; f)-D t -2t- (10f 4 + 2f 3 - 5t - 4), 
64f 4 -/f ) (3;0 = (t - 1) • J (4) (3; t) ■ D t 

-It- (448 f 6 + 32 f 5 + 95 f 4 - 220 1 2 - 1 12 1 - 128), 
576 ■ f 5 • jf>(4; t) = (f - 1) ■ /< 4) (4; f) ■ D, 

- 2 f • (5184 f 8 + 192 f 7 + 406 ? 6 + 1 148 ? 5 - 2471 f 3 - 1288 f 2 - 864 1 - 1 152), 

16384 • f 6 ■ 7j 4) (5; f) = (f - 1) ■ ^ 4) (5; f) • D, 

-It- (180224 f 10 + 4096 1 9 + 7488 f 8 + 15168 1 1 + 41307 1 6 

_ 83454 t 4 - 441 12 f 3 - 29952 1 2 - 22528 t - 32768). (F.2) 

Finally, the simplest jf\N; t), namely (2; ?) reads: 

16 f 6 ■ 4 4) (2; = 8 (f - l) 2 ■ J (4) (2; t) ■ D 2 

- t ■ (t - 1) • (432f 4 + 80f 3 - 99f 2 - 240f - 208) ■ D, 

+ 3 (1040r 5 - 1 176r 4 - 233? 3 - 100r 2 + 168? + 256). (F.3) 

Appendix G. Exact results for the C„ s 

The L(0' s have a new feature not previously seen. The inhomogeneous terms on the ODE's 
for C {2 ?(N; t) and C, (3) (N; f) begin at t N+a where a is 0,1 or 2 depending on the values of m. 
Therefore to the order needed for the polynomial solution the logarithms in the solution U2(N) 
never can contribute. However, for f^ N (N; t) the order of the inhomogeneous terms grows as 
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t 2N instead of t N . Therefore, since logarithms occur in u 2 (N) at order t N+l in order to find the 
polynomial solution to the 20-th order inhomogeneous equation in terms of the solutions of 
u 2 (N) and u x (N) we need to find linear combinations of solutions of the terms in the direct 
sum decomposition which cancel these logarithms. 

This procedure for solving the inhomogeneous equations is too cumbersome by itself to 
obtain explicit results as was done for f^ N (t) and f^ N (t). However, when the cancellation 
of logarithms is combined with the Wronskian cancellation method of section 5 it is 
possible to conjecture results for Cm\N; t) which have been verified to satisfy the 20th order 
inhomogeneous equations through N =10: 



(4) 2N+1 - (4 ) U 2 ( N + 1) 
—A, 



k; 



(4) 



c, 



(2) 



u 2 (N + 1) 



2 

+ 3' 



r 



(2) u\(N+ 1) 



o 



ul(N+l) „ +2 



o N -o 
- 2f3 N ■ Cg 

• f w+2 • F N+ i, 



( 2) u 2 2 (N + I) ■ u 2 (N) (2) u\{N+\) 

"5 -PN+l ' "5 



(G.l) 



u 3 JN + \)u 2 {N) 



o N 



2/3 N • q 



(2) u 2 (N + 1) • u 2 (N) , 2) u\{N+\) 



2 
+ - 

3 



u 2 (n + i)-u 2 2 (N) r oW2) « 3 (jV+l) 



+2C" 



2C 



Ni 2 



M 3 (N + 1) 



,2, »|(^+ 1) 

f 3 

(G.2) 



iV+1 



(4) 2W+1 jp(4) fl : 
C 2 - _6A: 'A 



2 "2 



M 2 (N + 1) • u 2 JN) 



^(2).,.2^ n ^(2) "2(A f + 1) ■ u 2 (N) | ^ + 1) 



2 
+ - 

3 



6/3 3 N ■ Cf ■ u\{N) ■ F N+l - 9f3 N ■ C\ 



(2) u 2 (N + 1) ■ u 2 (N) 



+ 6C 



(2) "2 



u 3 JN + 1) 



,2 "2 



M 2/ jV + 1) • u 2 (N) 



■t N+l -F N +3-P] 



u 2 (N + 1) • u 2 (N) 



2 y "' JV+2 P 

(G.3) 



r (4) 2i V_ + 2 ?(4 ). fi3 "2(A^ + 1) ; U 2 (N) 



2 
+ 3 



1 Pn'^q ■U 2 (Jy)-r N -2/Jjy-Cj ■ u 2 (iM ) • r N+ i -op N -L 2 ■ 3 /-at 



A^ 2 



(G.4) 
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(4 ) 2N+3 w 4 ) m a, „ ? (4) - 1 fl2 r (2) 2 rAA 





■F N +P 2 N -Cf -u 2 2 (N) 



(G.5) 



In order to construct the full C m , the expressions above are series expanded up to the order of 
palindromy, with palindromy determining the rest of the terms. The palindromy points of the 

are given as follows: m = : 2N + 1, m = 1 : 2N + 1, m = 2 : 2N + 2, m = 3 : 2N + 2, 
m — 4 : 2N + 3. Therefore, the expressions above give all terms to all except for the 
middle term of at order 2N + 2, which is determined such that all terms in cancel up 
to and including 2N + 3. 

Note that while these C$ guarantee that all terms will vanish up to and including 2N + 3, 
it is not obvious that the expansion at order 2N + 4 will match the expansion of even 
though it is the case. 
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